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One of the key challenges in spectroscopy is inhomogeneous broadening that masks 
the homogeneous spectral lineshape and the underlying coherent dynamics. A variety 
of techniques including four-wave mixing and spectral hole-burning are used in optical 
spectroscopj^ffl^ while in nuclear magnetic resonance (NMR) spin-echcP is the most 
common way to counteract inhomogeneity. However, the high-power pulses used in 
spin-echo and other sequenced®^] often create spurious dynamic^^ obscuring the sub¬ 
tle spin correlations that play a crucial role in quantum information applicationdSEEHIIl. 
Here we develop NMR techniques that allow the correlation times of the fluctuations 
in a nuclear spin bath of individual quantum dots to be probed. This is achieved 
with the use of frequency comb excitation which allows the homogeneous NMR line- 
shapes to be measured avoiding high-power pulses. We And nuclear spin correlation 
times exceeding 1 s in self-assembled InGaAs quantum dots - four orders of magnitude 
longer than in strain-free HI-V semiconductors. The observed freezing of the nuclear 
spin fluctuations opens the way for the design of quantum dot spin qubits with a 
well-understood, highly stable nuclear spin bath. 

Pulsed magnetic resonance is a diverse toolkit with applications in chemistry, biology and 
physics. In quantum information applications, solid state spin qubits are of great interest and are 
often described by the so called central spin model, where the qubit (central spin) is coupled to 
a fluctuating spin bath (typically interacting nuclear spins). Here microwave and radio-frequency 
(rf) magnetic resonance pulses are used for the initialization and readout of a qubilpSl, dynamic 
decoupling and dynamic controP of the spin bath. 

However, the most important parameter controlling the central spin coherenc^^^^^^ - the cor¬ 
relation time Tc of the spin bath fluctuations is very difficult to measure directly. The Tc is deter¬ 
mined by the spin exchange (flip-flops) of the interacting nuclear bath spins. By contrast pulsed 
NMR reveals the spin bath coherence time T 2 , which characterizes the dynamics of the transverse 
nuclear magnetizationPEIHl and is much shorter than Tc- The problem is further exacerbated in 
self-assembled quantum dots where quadrupolar effects lead to inhomogeneous NMR broadening 
exceeding 10 MHz (Refs.^^^^, making rf field amplitudes required for pulsed NMR practically 
unattainable. 

Here we develop an alternative approach to NMR spectroscopy: we measure non-coherent 
depolarisation of nuclear spins under weak noise-like rf fields. Contrary to intuitive expectation, 
we show that such measurement can reveal the full homogeneous NMR lineshape describing the 
coherent spin dynamics. This is achieved when rf excitation has a frequency comb profile (widely 
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FIG. 1. Prequency-comb technique for homogeneous NMR lineshape measurement, a, Timing 
diagram of the experimental cycle consisting of nuclear spin optical pumping (duration tpump), frequency- 
comb rf excitation (trf), and optical probing of the the nuclear spin state (tprobe)- b, The green line shows 
schematically a quantum dot NMR spectrum consisting of a central transition (CT) peak and two satellite 
transition (ST) bands. The inhomogeneous lineshape (width Az/inh) is a sum of a large number (>1000) 
of nuclear spin transitions with homogeneous linewidths Azxhom (shown with red lines). The thin black 
line shows the spectrum of the frequency-comb excitation with comb period /cp and total width Ar/comb 
exceeding Az/jnh. c, and d, demonstrate how experiments with varying /cp can reveal the width Az^hom of 
the NMR homogeneous lineshape (red lines). When Az/hom > /cp (c) all individual nuclei are uniformly 
excited by the frequency comb (shown with black lines). In the opposite case Az^hom < /cp (d) some of 
the nuclei (dashed line) are not excited resulting in a slow-down of nuclear spin dynamics. The transition 
between the two cases takes place when fcp Az/hom) allowing Az^ijojjj to be measured. 


used in precision optical metrologj^^SI) , \Ye then exploit non-resonant nuclear-nuclear interactions: 
the homogeneous NMR lineshape of one isotope measured with frequency comb NMR is used as a 
sensitive non-invasive probe of the correlation times Tc of the nuclear flip-flops of the other isotope. 
While initial studie^SEIISl suggested Tq ~ 100 fis for nuclear spins in III-V semiconductors, it was 
recently recognizecP^^^ that quadrupolar effects may have a significant impact in self-assembled 
quantum dots. Here we for the first time obtain a quantitative measurement of extremely long 
Tc > 1 s revealing strong freezing of the nuclear spin bath - a crucial advantage for quantum 
information applications of self-assembled quantum dots. 

The experiments were performed on individual neutral self-assembled InGaAs/GaAs quantum 
dots at magnetic field Hz = 8 T. All measurements of the nuclear spin depolarisation dynamics 
employ the pump-depolarise-probe protocol shown in Fig. dh- Here we exploit the hyperfine 
interaction of the nuclei with the optically excited electroi]^^^^^ both to polarise the nuclei (pump 
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pulse) and to measure the nuclear spin polarisation in terms of the Overhauser shift AE'hf in the 
QD photoluminescence spectrum (probe pulse). The rf magnetic field depolarising nuclear spins 
is induced by a small copper coil. (Further experimental details can be found in Methods and 
Supplementary Note 1.) 

All isotopes in the studied dots possess non-zero quadrupolar moments. Here we focus on the 
spin I = 3/2 nuclei ^^Ga and ^^As. The strain-induced quadrupolar shifts result in an inhomoge- 
neously broadened NMR spectruirP^II a,s shown schematically by the green line in Fig. [I}3. The 
spectrum consists of a central transition (CT) —1/2 -H- -1-1/2 and two satellite transition (ST) 
±1/2 -H- ±3/2 peaks. The NMR spectrum with inhomogeneous linewidth Az/inh consists of indi¬ 
vidual nuclear spin transitions (shown with red lines) with much smaller homogeneous linewidth 
Az^hom- 

To make a non-coherent depolarisation experiment sensitive to the homogeneous NMR lineshape 
rf excitation with a frequency comb spectral profile is used. As shown in Fig. (black line) the 
frequency comb has a period of /cp and a total comb width Ai/comb exceeding Az^inh- The key 
idea of the frequency comb technique is described in Figs. [TJ: and d where two possible cases are 
shown. If the comb period is small (/cp < Az'hom) Fig- all nuclear transitions are excited by 
a large number of rf modes. As a result all nuclear spins are depolarised at the same rate and we 
expect an exponential decay of the total nuclear spin polarisation. In the opposite case of large 
comb period (/cp > Aphom) Fig. [^) some of the nuclear transitions are out of resonance and are 
not excited (e.g. the one shown by the dashed red line). As a result we expect a slowed-down 
non-exponential nuclear depolarisation. The experiments are performed at different /cp; the /cp 
for which a slow-down in depolarisation is observed gives a measure of the homogeneous linewidth 

^hom* 

Experimental demonstration of this technique is shown in Fig. [2jr. The Overhauser shift 
variation AFi^f of ^^Ga is shown as a function of the depolarising rf pulse duration tj-f for different 
/cp. For small /cp =80 and 435 Hz an exponential depolarisation is observed. However, when 
/cp is increased the depolarisation becomes non-exponential and slows down dramatically. The 
detailed dependence /cp) is shown as a colour-coded plot in Fig. Ep- The threshold value 

of /cp (marked with a white arrow) above which the nuclear spin dynamics becomes sensitive 
to the discrete structure of the frequency comb, provides an estimate of Aphom ~ 450 Hz. Such 
a small homogeneous linewidth is detected in NMR resonances with inhomogeneous broadening 
of Al/inh r\j 6 MHz (Ref.^ demonstrating the resolution power of frequency-comb non-coherent 
spectroscopy. 


4 




b 



rf pulse duration, (s) 


FIG. 2. Measurement of the homogeneous NMR lineshape in self-assembled quantum dots 
using frequency-comb excitation, a, The change in the polarisation of the nuclear spins (measured 
in terms of the change in the Overhauser shift Aifhf) is shown as a function of the rf pulse duration Gf 
(symbols) at = 8 T and different comb periods fcp- Lines show model fitting (see text), b, A full 2D 
plot of Ai?hf as a function of ipp and /cp in the same experiment as in a. There is a clear slow-down of 
the nuclear spin depolarisation at /cp > 450 Hz (shown with a white arrow) providing an estimate of the 
homogeneous linewidth of Azzpom- 


The information revealed by frequency-comb spectroscopy is not limited to linewidth estimates. 
An accurate determination of the full homogeneous lineshape is achieved with modeling based 
on solving an integral equation (see details in Methods and Supplementary Note 2). We use the 
following two-parameter phenomenological model for the homogeneous lineshape: 

L{u) (X (^1 + 4:{V2-1 )-^—^ , (1) 

where Azz^om is the homogeneous full width at half maximum and A: is a roll-off parameter that 
controls the tails of the lineshape (the behavior of L{u) at large u). For k = 1 the lineshape 
corresponds to Lorentzian, while for A: —)• oo it tends to Gaussian; in this way Eq. seamlessly 
describes the two most common lineshapes. Using Az/hom and k as parameters we calculate the 
model AEhf(trf, /cp) dependence and fit it to the experimental AEhf(trf) /cp) to find an accurate 
phenomenological description of the homogeneous NMR lineshape in self-assembled quantum dots. 

The solid line in Fig. shows the best-fit lineshape (Ai/hom ~ 221 Hz and k ~ 1.67) for 
the measurement shown in Figs. Ifr, b. The dashed and dashed-dotted lines in Fig. show 
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FIG. 3. Homogeneous lineshape modeling, a, b, Model homogeneous NMR lineshapes of nuclei 
shown on linear (a) and logarithmic (b) scale. Solid lines show the best-fit lineshape with a full width at 
half maximum A^'hom ~ 221 Hz and a roll-off parameter k « 1.67. Dashed and dashed-dotted lines show for 
comparison Lorentzian and Gaussian lineshapes with the same A^'hom- c-e, The calculated Ai?hf(trf, /cp) 
dependencies for the lineshapes in (a, b). For the fitted lineshape (c) an excellent agreement with the 
experiment in Fig. is found, while calculations with Lorentzian (d) and Gaussian (e) lineshapes give 
markedly different results demonstrating the sensitivity of the frequency comb technique. 

for comparison the Lorentzian (k = 1) and Gaussian {k —)> oo) lineshapes with the same Az/hom- 
The difference in the lineshape tails is seen clearly in Fig. where a logarithmic scale is used. 
The model A£'hf(trf)/cp) dependence calculated with the best fit parameters is shown in Fig. 
and with lines in Fig. [2}r - there is excellent agreement with experiment. By contrast modelling 
AFihf(^rf 5 /cp) with Lorentzian (Fig. [^) and Gaussian (Fig. [^) lineshapes show a pronounced 
deviation from the experiment, demonstrating the excellent sensitivity of the frequency-comb spec¬ 
troscopy to accurately probe the homogeneous spectral lineshape. 

We have also performed frequency comb NMR spectroscopy on ^^As nuclei (Fig. |^). Despite 
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their larger inhomogeneous broadening Az^inh ~ 18 MHz the model htting reveals even smaller 
At'hom ~ 117 Hz and k ~ 1.78. The frequency-comb measurements are in agreement with the 
previous findings based on spin-echo NMR measurement^ indeed, from Ai/hom derived here we 
can estimate the nuclear spin coherence time T 2 ~ l/(7rAr'hom) ~ 1-4 and 2.7 ms for ^^Ga and ^®As, 
in good agreement with the corresponding spin-echo T 2 ~ 1.2 and 4.3 ms. On the other hand, spin- 
echo could only be measured on central transitions for which Ai/inh is relatively small. Moreover 
pulsed NMR does not allow determination of the full homogeneous lineshape, which for dipole- 
dipole interactions typically has a ’’top-hat”-like (Guassian) profil^^fi]^ And most importantly, due 
to the parasitic effects such as ” instantaneous diffusion”® and spin locking pulsed NMR does not 
reveal the characteristic correlation time Tc of the spin exchange (spin flip-flop) between the nuclei 
in the absence of rf excitation. 

As we now show, the non-Gaussian lineshapes can be understood and can be derived using 
experiments with two frequency combs exciting nuclei of two isotopes (^^As and ^^Ga). The two- 
comb experiment is similar to that shown in Fig. ir: we excite ’^^As nuclei with a frequency comb 
to measure their homogeneous lineshape. The difference is that now we simultaneously apply a 
second comb exciting the ^^Ga spins. Importantly, in this experiment the ^^Ga nuclei are first 
fully depolarised after the optical nuclear spin pumping -- in this way the excitation of ^^Ga has 
no direct effect on the measured hyperfine shift AHhf. By contrast it leads to ’’heating” of the 
^^Ga spins which has only an indirect effect on by changing the ^^As lineshape via dipolar 

coupling between ^^Ga and ^®As spins. The result of the two-comb experiment is shown in Fig. 

a clear increase of Ai/hom for ^®As is observed. From model fitting we find that ^^Ga ’’heating” 
leads to a 3 times broader homogeneous linewidth Ai^hom ~ 355 Hz of ^^As and its homogeneous 
lineshape is modified towards Gaussian, observed as increased k ^ 2.32. 

To explain this result we note that the NMR lineshape is a statistical distribution of NMR 
frequency shifts of each nucleus produced by its dipolar interaction with all possible configurations 
of the neighboring nuclear spins. However, the frequency comb experiment is limited in time (up 
to ~100 s as shown in Figs. |^, b). If the nuclear spin environment of each ^®As nucleus does 
not go through all possible configurations during the measurement time, the frequency shifts are 
effectively static, and hence are eliminated from the lineshape as for any other inhomogeneous 
broadening. 

Thus we conclude that the narrowed, non-Gaussian (k ~I.6-1.8) homogeneous NMR lineshape 
arises from the ’’snapshot” nature of the frequency comb measurement, probing the strongly frozen 
nuclear spin configuration. When the additional ^^Ga ’’heating” excitation is applied it ’’thaws” 
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FIG. 4. Probing the correlation times of the nuclear spin bath fluctuation, a, b, Frequency-comb 
measurement on nuclei at B^, =8 T without (a) and with (b) additional frequency-comb excitation 
’’heating” the ^^Ga nuclei. The increase of A^'hom of ^^As under the ^^Ga ’’heating” reveals the strongly 
suppressed nuclear spin fluctuations of ^^Ga. c, rf-induced depolarisation time tas of ^®As nuclei at fixed 
/cp = 1.47 kHz as a function of the amplitude /3i of the addition rf excitation ’’heating” either the central 
transition (CT, triangles) or the satellite transition (ST, squares) of ^^Ga. The bottom scale shows /?i 
expressed in terms of the rf-induced spin-flip time rca of ^^Ga. The rca at which a marked decrease of tas 
is observed corresponds to the correlation time Tc of the ^^Ga nuclear spin fluctuations: By extrapolating 
the power-law regions (dashed lines) to the intersection with the maximum tas (solid line) we deduce 
"^c.CT ~ 1 sec and Tc^st ~ 20 sec for the CT and ST transitions. 
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the spins, detected as broadening of the ’^^As lineshape (as demonstrated in Figs. [^, b). We 
use such sensitivity of the ^®As lineshape to measure the dynamics of the ^^Ga equilibrium spin 
bath fluctuations. Based on the results of Fig. |^the ^^As spins are now excited with a frequency 
comb with a fixed /cp = 1.47 kHz for which the ^®As depolarisation dynamics is most sensitive to 
the ^^Ga ’’heating”. Furthermore, we now use selective ’’heating” of either the —1/2 -f-)- +1/2 CT 
or the +1/2 •(-)> +3/2 ST of ^^Ga. The amplitude /3i of the ’’heating” frequency comb is varied - 
the resulting dependencies of the ^^As depolarisation time tas are shown in Fig. by the squares 
and triangles for CT and ST ’’heating” respectively. For analysis we also express /3i in terms of 
the rf-induced ^^Ga spin-flip time rca (bottom scale in Fig. [^, see details in Methods). 

It can be seen that for vanishing "^^Ga excitation (/3 —)> 0) the '^^As depolarisation time tas 
is constant. In this ’’frozen” regime the rf-induced spin-flip time rca of ^^Ga is larger than the 
correlation time Tq of the ^^Ga intrinsic spin flip-flops (rca > +;)• As a result tas is determined 
only by the rf excitation of ^®As itself. However, when (3i is increased to ~1- 10 nT Hz the 
rf induced spin-flips of ^^Ga nuclei become faster than their intrinsic flip-flops (rca < +)• Such 
’’thawing” of ^^Ga broadens the ^^As lineshape (via heteronuclear interaction), and is observed as 
a reduction of tas- Thus the transition from the ’’frozen” to ’’thawed” regimes takes place when 
'T'Ga ~ +) allowing Tc to be determined. In Fig. we extrapolate graphically (dashed lines) the 
power-law dependence in the ’’thawed” regime. The points of the intersections with the limiting 
value of tas in the ’’frozen” regime (solid horizontal line) yield correlation times Tc^ct ~ 1 s for CT 
and Tc,sT ~ 20 s for ST. 

The observed + > 1 s exceeds very strongly typical nuclear dipolar flip-flop times in strain-free 
HI-V solids Tc ~ 100 We attribute the extremely long Tc in self-assembled quantum dots 

to the effect of inhomogeneous nuclear quadrupolar shifts making nuclear spin flip-flops energeti¬ 
cally forbiddenP^II, This interpretation is corroborated by the observation of Tc^st +,CT, since 
quadrupolar broadening of the ST transitions is much larger than that of the CT^l^. Furthermore, 
the ^^Ga spins examined here have the largest gyromagnetic ratio 7 and the smallest quadrupolar 
moment Q, so we expect that all other isotopes in InGaAs have even longer Tc, resulting in the 
overall Tc > 1 s of the entire quantum dot nuclear spin bath. This implies that in high magnetic 
fields the spin-echo coherence times of the electron and hole spin qubits in self-assembled dots are 
not limited by the nuclear spin bath up to sub-second regime^^^^oIiSIlZl Provided that other mech¬ 
anisms of central spin dephasing, such as charge fluctuationd^sHUl ^re eliminated, this would open 
the way for optically active spin qubit networks in HI-V semiconductors with coherence properties 
previously achievable only in nuclear-spin-free materialJ^^®. 
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Since the frequency-comb technique is not limited by artifacts in the spin dynamics hampering 
pulsed magnetic resonance, it allows detection of very slow spin bath fluctuations. Such sensi¬ 
tivity of the method can be used for example to investigate directly the effect of the electron or 
hole on the spin bath fluctuations in charged quantum dots, arising for example from hyperfine- 
mediated nuclear spin interactions. The experiments can be well understood within a classical rate 
equation model, while further advances in frequency comb spectroscopy can be expected with the 
development of a full quantum mechanical model. Furthermore, the simple and powerful ideas of 
frequency-comb NMR spectroscopy can be readily extended beyond quantum dots: as we show in 
Supplementary Note 3 the only essential requirement is that the longitudinal relaxation time Ti 
should be larger (by about two orders of magnitude) than the transverse relaxation time T 2 , which 
is usually the case in solid state spin systems. Finally our approaches in the use of frequency combs 
can go beyond NMR, and for example enrich the techniques in optical spectroscopy. 


I. METHODS SUMMARY 

Sample structures and experimental techniques The experiments were performed on 
individual neutral self-assembled InGaAs/GaAs quantum dots. The sample was mounted in a 
helium-bath cryostat (r=4.2 K) with a magnetic field Rz = 8 T applied in the Faraday configura¬ 
tion (along the sample growth and light propagation direction Oz). Radio-frequency (rf) magnetic 
field Rrf perpendicular to Rz was induced by a miniature copper coil. Optical excitation was used 
to induce nuclear spin magnetization exceeding 50%, as well as to probe it by measuring hyperfine 
shifts in photoluminescence spectroscopj®. 

Two sample structures have been studied, both containing a single layer of InGaAs/GaAs 
quantum dots embedded in a weak planar microcavity with a Q-factor of ~250. In one of the 
samples the dots emitting at ~ 945 nm were placed m a p — i — n structure, where application of a 
large reverse bias during the rf excitation ensured the neutral state of the dots. The results for this 
sample are shown in Fig. The second sample was a gate-free structure, where most of the dots 
emitting at ~ 914 nm are found in a neutral state, although the charging can not be controlled. 
Excellent agreement between the lineshapes of both ^^Ga and ^®As in the two structures was found, 
confirming the reproducibility of the frequency-comb technique. 

Homogeneous lineshape theoretical model. Let us consider an ensemble of spin 1 = 1/2 
nuclei with gyromagnetic ratio 7 and inhomogeneously broadened distribution of nuclear resonant 
frequencies Unuc- We assume that each nucleus has a homogeneous absorption lineshape 
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with normalization L{v)dv = 1. A small amplitude (non-saturating) rf field will result in 
depolarisation, which can be described by a differential equation for population probabilities p±i /2 
of the nuclear spin levels = zbl/2 


d{p+ii 2 -P-il 2 )ldt = -bh(p+i /2 - P- 1 / 2 )• 


( 2 ) 


For frequency-comb excitation the decay rate is the sum of the decay rates caused by each rf mode 
with magnetic field amplitude Bi, and can be written as: 


VF(z^nuc) = X] - 1 ^ 1 - j7cp)/cP, 


2/cp 


(3) 


i=o 


where the summation goes over all modes with frequencies -l-j/cp is the frequency of 

the first spectral mode). 

The change in the Overhauser shift produced by each nucleus is proportional to p+i/ 2 ~P-i /2 
and according to Eq. has an exponential time dependence oc exp(—fF(t'nuc)i)- The quantum 
dot contains a large number of nuclear spins with randomly distributed absorption frequencies. 
Therefore to obtain the dynamics of the total Overhauser shift we need to average over t'nuo which 
can be done over one period /cp since the spectrum of the rf excitation is periodic. Furthermore, 
since the total width of the rf frequency comb Ai/comb is much larger than /cp and Aphom; the 
summation in Eq. [^can be extended to Too. Thus, the following expression is obtained for the time 
dependence AEhf(t,/cp), describing the dynamics of the rf-induced nuclear spin depolarisation: 


AEhf(t, /cp) ^1^-1 
AEhf(t —)• oo) 


J exp I X] - jfcp)fcP I rf^'nuc- 

0 \ i=-oo / 


(4) 


Equation 1^ describes the dependence AEhf(t,/cp) directly measurable in experiments such as 
shown in Fig. A£'iif(t —>■ oo) is the total optically induced Overhauser shift of the studied 
isotope and is also measurable, while fcp and Bi are parameters that are controlled in the ex¬ 
periment. We note that in the limit of small comb period fcp —)• 0 the infinite sum in Eq. 
tends to the integral L{v)di' = 1 and the Overhauser shift decay is exponential (as observed 
experimentally) with a characteristic time 


r = 2fcp/h^Bf) (5) 

Equation is a Fredholm’s integral equation of the first kind on the homogeneous lineshape 
function L{v). This is an ill-conditioned problem: as a result finding the lineshape requires some 
constraints to be placed on L{v). Our approach is to use a model lineshape of Eq. After 
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substituting L{u) from Eq. the right-hand side of Eq. [^becomes a function of the parameters 
At'hom and k which we then find by least-squares fitting of Eq. |^to the experimental dependence 
AEhf(t, /cp). 

This model is readily extended to the case of / > 1/2 nuclei. Eq. [^becomes a tri-diagonal system 
of differential equations, and the solution (Eq. contains a sum of multiple exponents under the 
integral. These modifications are straightforward but tedious and can be found in Supplementary 
Note 2. 

Derivation of the nuclear spin bath correlation times. Accurate lineshape modeling is 
crucial in revealing the ^^As homogeneous broadening arising from ’^^Ga ’’heating” excitation (as 
demonstrated in Figs. |^, b). However, since a measurement of the full A£'hf(t,/cp) dependence 
is time consuming, the experiments with variable ^^Ga excitation amplitude f3i (Fig. i:) were 
conducted at fixed /cp = 1.47 kHz exceeding noticeably the ^^As homogeneous linewidth Aphom ~ 
117 Hz. To extract the arsenic depolarisation time tas we fit the arsenic depolarisation dynamics 
AEhf(t) with the following formulae: AEhf(trf) = AEhf(trf —^ oo)(l — exp[—(Gf/'TAs)’'])) using r as 
a common fitting parameter and tas independent for measurements with different roa- We find 
r ~ 0.57, while the dependence tas on rca obtained from the ht is shown in Fig. with error 
bars corresponding to 95% confidence intervals. 

The period of the ^^Ga ’’heating” frequency comb is kept at a small value /cp = 150 Hz 
ensuring uniform excitation of all nuclear spin transitions. The amplitude of the ’’heating” comb 
is defined as //i = where Bi is magnetic field amplitude of each mode in the comb 

(further details can be found in Supplementary Note 1). To determine the correlation times we 
express /3i in terms of the rf-induced spin-flip time rca- The rca is defined as the exponential 
time of the ^^Ga depolarisation induced by the ’’heating” comb and is derived from an additional 
calibration measurement. The values of /3i shown in Fig. correspond to the experiment on the 
CT and are calculated using Eq. [^as where 7 is the ^^Ga gyromagnetic ratio and 

TQa is experimentally measured. The additional factor of 4 in the denominator is due to the matrix 
element of the GT of spin I = 3/2. For experiments on ST the /3i values shown in Fig. must be 
multiplied by i/4/3. 
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Supplementary Figure 1. Frequency combs for NMR experiments, a, Schematic NMR spectrum of 
a strained self-assembled InGaAs quantum dot at = 8 T. Sharp peaks arise from the central transitions 
(CTs) between —1/2 and -Fl/2 nuclear spin levels of each isotope. There are two satellite transitions (STs) 
for the spin-3/2 isotopes of Ga and As, and eight STs for the spin-9/2 indium. Inhomogeneous quadrupolar 
broadening results in significant overlap of the ^^^In and ®®Ga inhomogeneous lineshapes, while the overlap 
of ^®As and ^^^In is minimal. The ^^Ga resonance is well isolated, b-d, Schematic spectra of the frequency 
combs used in experiments. The ^^As-full comb is used to excite the entire inhomogeneous lineshape of ^®As 
with minimum effect on Similarly ^^Ga-full comb excites the entire ^^Ga resonance. Frequency comb 

^^Ga-ST is used to excite selectively the -Fl/2 e-i- -F3/2 ST of ^^Ga, while the comb ’’^Ga-GT is used to excite 
the —1/2 -44 -F1/2 GT (zoom in is shown in d). The parameters of all frequency combs are summarized in 
[Supplementary Table l[ further details can be found in [Supplementary Note 1 A| 
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Supplementary Figure 2. Time diagram of an optically detected frequency-comb nuclear mag¬ 
netic resonance experiment. See detailed explanation in [Supplementary Note 1^ 



Supplementary Figure 3. Experimental verification of the rate equation model, a, Dynamics of the 
^^Ga nuclear spin depolarisation at = 8 T induced by frequency comb rf excitation with comb period 
/cp = 31.25 Hz. The results at high rf field density j3i = 86.0 nT Hz“^/^ (low rf field density /3i = 24.2 nT 
are shown by the squares (triangles). Solid lines show exponential fits with decay time r « 23 ms 
(r « 290 ms) for high (low) rf field density, b, c, Full AEhf(t, /cp) dependence measured at high (b) and 
low rf (c) field densities. 
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TABLE Supplementary Table 1. Frequency comb parameters used in NMR experiments. 


a, Homogeneous lineshape experiments (Fig. 2): 

Isotope 

'^®As 

Ga 



Frequency comb 

75As-full 

7iGa-full 



Comb central frequency t'comb.o (MHz) 

58.81 

104.80 



Comb spectral width At'comb (MHz) 

18 

9 



Comb period fcp (Hz) 

varied 

varied 



Rf field density f3i (nT Hz“^/^) 

65.5 

39.1 



b, Lineshape broadening experiments (Figs. 

4a, 4b): 




Isotope 

'^®As 

71 Ga 



Frequency comb 

^^As-full 

7iGa-full 



Comb central frequency :^comb,o (MHz) 

58.81 

104.80 



Comb spectral width Ar'comb (MHz) 

18 

8 



Comb period /cp (Hz) 

varied 

159 



Rf field density Pi (nT Hz“^/^) 

30.5 

27.6 



c, Correlation time experiments (Fig. 4c): 

Isotope 

'^®As 

71 Ga 

7iGa 

7iGa 

Frequency comb 

^^As-full 

7iGa-full 

7iGa-ST 

7iGa-CT 

Comb central frequency t'comb.o (MHz) 

58.81 

104.80 

106.10 

104.80 

Comb spectral width At'comb (MHz) 

18 

8 

2.5 

0.05 

Comb period fcp (Hz) 

1466 

159 

150 

150 

Rf field density /?i (nT Hz“^/^) 

30.5 

27.6 

varied 

varied 
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Supplementary Note 1. DETAILS OF EXPERIMENTAL TECHNIQUES 


A. Frequency combs 

The key novel findings of this work are based on the use of radiofrequency (rf) excitation with 
a frequency comb spectral profile. Here we give detailed parameters of the frequency combs used 
in the experiments on self-assembled quantum dots. 

Supplementary Figure shows a schematic NMR spectrum of a self-assembled InGaAs quan¬ 
tum dot at = 8 T, based on results of the inverse NMR measurementFour sharp peaks 
arise from the central transitions (CTs) —1/2 -H- -1-1/2 of each isotope. All spin-3/2 nuclei (^®As, 
®®Ga and ’^^Ga) have two satellite transitions (STs) —3/2 •(-)> — 1/2 and -1-1/2 -H- -1-3/2 observed 
as inhomogeneously broadened bands on both sides of the GTs. The spin-9/2 ^^®In nuclei have 
a total of eight STs. For clarity these are shown as four bands on each side of the CT, although 
in experimental NMR spectra the peaks arising from different STs merge forming two inhomoge¬ 
neously broadened bands on both sides of the ^^®In CT^^. As shown in Supplementary Figure 
the spectral contributions from the ^^®In and ®®Ga nuclei overlap significantly; furthermore, there 
is a small overlap between the ^®As and ^^^In NMR resonances. 

In the experiments we use rf frequency combs that are designed to influence only the chosen 
transition (s) of one isotope without affecting the other isotopes as demonstrated in Supplementary 
Figures [Ija-d. The frequency-comb labeled ^^Ga-full has a total width At'comb of 8 or 9 MHz. This 
comb entirely covers the inhomogeneous lineshape of ^^Ga and thus uniformly excites all nuclear 
spin transitions of this isotope. Since ^^Ga has a large resonance frequency the ^^Ga-full comb 
does not affect the polarisation of the other isotopes. To achieve selective excitation of the ^^Ga 
ST (-1-1/2 <-)■ -1-3/2) we use frequency comb ^^Ga-ST with Ai/comb = 2.5 MHz (Supplementary Fig. 
[^). Similarly, selective excitation of ^^Ga GT is achieved with a narrow comb ^^Ga-GT with a 
width of Ai/comb = 50 kHz as shown in Supplementary Figs. P, d. 

In the case of ^®As we use a frequency comb ^^As-full that excites the entire inhomogeneouse 
resonance line. This comb inevitably excites some of the ^^®In nuclear transitions, mostly —9/2 -H- 
—7/2 STs. However, excitation of the —9/2 -f-)- —7/2 ST alone has a negligible effect on the overall 
change in the nuclear polarisation of the spin-9/2 nuclei of ^^^In. Thus the ^®As-full comb with 
the optimum width Ai/comb = 18 MHz is used in experiments for selective depolarisation of ^^As. 

The central frequencies r'comb, 0 ) the widths At'comb and the comb periods (spectral separation 
between the adjacent modes) /cp of all the combs used in experiments are listed in [Supplementary 
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Table 1 Depending on the experiment the fcp is either varied or kept constant. Typically the 


values of /cp ranging from 30 Hz to 21 kHz are employed, so that the total number of modes in 
the comb = Az/comb/ZcP + 1 ranges from ~330 to ~600000. 


The phases of individual modes of the frequency comb are chosen in a way that minimizes 
the peak power for a given average power (i.e. a waveform with the minimum crest-factor). In 
particular the following expression satishes this criterion: 


Nm 

Brf{t) = Bi COS 


{ui + {j - l)/cp) t + TT 


iVm y’ 


( 1 ) 


where the summation goes over all modes, and vi is the frequency of the first mode of the comb. In 
experiments the frequency-comb signal is generated by an arbitrary waveform generator equipped 
with a 64 million points memory. 

Each mode of the frequency comb has the same amplitude Bi of the rf oscillating magnetic field. 
In experiments where the comb period /cp is varied Bi has to be adjusted to maintain the same 
total power of the frequency comb excitation. Since the power is proportional to B^ the H^//cp 
ratio has to be kept constant. Thus the amplitude of the frequency comb can be conveniently 
characterized by the magnetic field density j3i = Bi /VTcP • 

The amplitude of the frequency comb rf magnetic field can be calibrated from an additional 
pulsed Rabi oscillation experiment on the CT of a selected isotopd^. The Rabi oscillation circular 
frequency WRabi can be measured experimentally and is proportional to the rf held amplitude Rj-f 
in the rotating frame: 


Rrf = 


‘^Rabi 
27 ’ 


( 2 ) 


where 7 is nuclear gyromagnetic ratio and the factor 2 originates from the dipole transition matrix 
element of the CT of a spin-3/2 nucleud^. Throughout all experiments, we monitored the applied 
rf helds via a pick-up coil that was placed close to the sample and connected to a spectrum analyzer. 
By comparing the voltages induced by the frequency comb modes with the voltage associated with 
a held Rrf in the Rabi oscillation experiment, we derived the values of Bi and /3i of the comb. The 


values of (3i used in different experiments are given in Supplementary Table 1 - these correspond 


to the rotating frame, i.e. the physical values of the magnetic held induced by the coil are twice 
as large. 
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B. Optical pump-probe techniques for frequency-comb NMR 


Supplementary Fig. shows the time sequence of optical and rf excitation pulses used in a 
frequency comb measurement of the equilibrium nuclear spin bath fluctuations of the ^^Ga isotope 
(the results are shown in Fig. 4c of the main text). As explained in the main text this experiment 
is based on measuring the rf induced dynamics of the ^^As spins in the presence of additional rf 
excitation of ^^Ga spins. The experimental cycle consists of the following four stages described 
below. 

Optical nuclear spin pumping. At the start of each new measurement cycle, the nuclear spin 
bath is reinitialized optically. This is achieved with optically induced dynamic nuclear polarisation 
(DNP)IS3IS11 under high power, circularly polarised laser excitation, spin polarised electrons are 
created. These electrons can efficiently transfer their polarisation to the nuclear spin bath via 
hyperfine interactiorP^. We use a~ polarised excitation with a pump laser operating at ~850 nm, in 
resonance with the QD wetting layer. By using sufficiently long pumping times (tpump = 5.5 —6.5 s) 
and high powers ~ lOPgat (where Pgat is the saturation power of the QD ground states), we create 
a reproducible and high degree of nuclear polarisation. 

Depolarisation of Ga isotope. Optical spin pumping polarises the nuclei of all isotopes. How¬ 
ever, in order to probe the equilibrium fluctuations of ^^Ga spins their longitudinal relaxation has 
to be excluded from the measured dynamics. For that ^^Ga nuclear polarisation has to be erased, 
which is achieved by exciting the spins with a ^^Ga-full frequency comb for a sufficiently long time 
^Depoi = 1.2 s. To simplify experimental implementation the depolarising rf is kept on during the 
nuclear spin pumping stage as well, which has no effect on the experimental results. 

Frequency comb rf excitation. Following the spin bath preparation (optical DNP and ^^Ga 
depolarisation), the main rf excitation (variable duration trf) is applied. For the spin bath fluctu¬ 
ation measurement (Fig. 4c) this excitation is a sum of the ^®As-full frequency comb and either 
the ^^Ga-ST or the ^^Ga-CT comb. Since ^^Ga is completely depolarised by the previous pulse, 
all changes in the total nuclear polarisation at this stage are solely due to the ^®As depolarisation. 
In this way we ensure that it is the depolarisation dynamics of ^®As that is measured, while the 
^^Ga-ST or ^^Ga-GT ’’heating” excitation only induces nuclear spin-flips of the corresponding ^^Ga 
transition. 

Optical probing of the nuclear spin polarisation. At the end of the experiment cycle, a short 
probe laser pulse is applied and the resulting photoluminescence spectrum is collected by a 1 m dou¬ 
ble spectrometer with a GGD. The changes in the quantum dot Zeeman splitting (the Overhauser 


6 




shift £'hf) are used to probe the nuclear spin state. The probe laser is non-resonant (~ 850 nm), 
yet unlike the pump laser it is linearly polarised and the probe power and duration tprobe are 
chosen such that no noticeable DNP is induced and the final nuclear spin polarisation is measured 
accurately. Typical probe parameters used in the experiments were tprobe = 4 ms and ~ Psat/10 
for QDs in the p — i — n diode sample and tprobe = 60 ms and ~ Psat/50 in the gate-free structure. 
Depending on the QD photoluminescence intensity the experiment cycle was repeated 10 — 40 times 
to achieve the optimum signal-to-noise ratio. 

For the purpose of data analysis we are interested in measuring the rf-induced change in the 
nuclear spin polarisation (rf-induced change in the Overhauser shift APhf); for that we perform a 
control measurement where ^®As-full comb is off, and subtract the resulting QD Zeeman splitting 
from the Zeeman splittings obtained in the measurements with As excitation. In this way = 
0 corresponds to no nuclear spin depolarisation induced by the rf. 

The diagram of Supplementary Fig. also describes the other types of frequency comb NMR 
measurements presented in the main text with the following modifications: For the line broadening 
measurements shown in Figs. 4a, b the main excitation is a sum of the ^^As-full and ^^Ga-full 
frequency combs (the ^^Ga-full is off for the measurement in Fig. 4a and is on for Fig. 4b). The 
homogeneous lineshape measurement (Fig. 2) is performed without the additional rf depolarisation 
excitation, while for the main rf excitation the ^^Ga-full comb is used. 


Supplementary Note 2. THEORETICAL MODEL FOR NUCLEAR SPIN DYNAMICS 
UNDER FREQUENCY-COMB EXCITATION 

The Methods section of the main text describes the model for spin 1 = 1/2 nuclei. Here we 
consider a more general case of nuclear spins I > 1/2 with gyromagnetic ratio 7 . We consider the 
Overhauser shifts of nuclei of only one isotope, which is justified since the polarisation of other 
isotopes stays constant during the measurement and can be neglected. In an external magnetic field 
the nuclear spin state is split into 27-1-1 states with spin projections = —7, —7-|-1, ... 7 — 1, +1. 
In our classical rate equation model we assume that each nuclear spin has a probability pm to be 
found in a state with = m with normalization condition 

I 

Pm{t) = 1. (3) 

m=—I 


7 


We also assume that at t = 0 optical pumping initializes all nuclear spins into a Boltzman distri¬ 
bution 


Pmit = 0 ) = pm,o oc exp(^m), 


(4) 


so that nuclear spins can be characterized by a temperature Tnuc oc 1/^. When optical nuclear spin 
pumping is used, Tnuc is very small compared to the spin bath temperature T, so the equilibrium 
nuclear polarisation can be neglected. 

Application of the radiofrequency (rf) excitation leads to the changes in population probabilities. 
The rf excites only dipole-allowed transitions for which changes by ±1. In the experiment we use 
weak, non-saturating radiofrequency fields. Thus instead of the full Bloch equations for nuclear 
magnetization, the evolution of the population probabilities of the states with = m can be 
described with the following first-order differential equation (see further details in 
Note 3): 


Supplementary 


dpm/dt — (Wm—l,m T ^^m,m+l)Pm{t) T IWn—1(0 4“ ^^m^m+lPm+lit) i / < Ul < / 
dpi/dt = Wi-ij{-pi{t) +pi_i{t)) (5) 

dp_j/dt = W-i-i+i{-p-i{t) +p_j+i{t)). 

Here the first equation describes the states with —I < Iz < I. Its first term is due to nuclei with 
Iz, = m making transitions into the = m — 1 and m-|-l states, whereas the second and third terms 
describe the opposite case of nuclei transitioning into the Iz = m state from the Iz = m — 1 and 
Iz = m -|- 1 states respectively. The second and third equations correspond to the case of m = —I 
and m = +I respectively. Taken for all m, for which —I < m < I, the Supplementary Eq. [^yields 
a system of 21 -|- 1 first-order ordinary differential equations (ODEs) for 21 -|- 1 time-dependent 
variables Pm{t) with initial conditions given by Supplementary Eq. Due to the normalization 
condition of Supplementary Eq. only 21 variables and equations are independent. This system 
of ODEs has a tridiagonal matrix with coefficients determined by the rf induced transition rates 
Wm,m+i which satisfy a symmetry condition Wm,m+i = Wm+i,m- 

Similar to the case of I = 1/2, each transition rate Wm,m+i resulting from the frequency- 
comb excitation is a sum of transition rates caused by individual rf modes each having magnetic 
field amplitude Bi. We assume that each nuclear transition has the same broadening described 
by the homogeneous lineshape function L{v), with normalization L{i')di' = 1. Due to the 
inhomogeneous quadrupolar shifts the NMR transition frequency t'm.m-i-i is generally different for 
each pair of spin levels B = m and Iz = m -|- 1 (even for one nucleus). Thus for the transition rates 




we can write: 


W, 


m,m+l 




m+l) — 


(/ — m)(I + m + 1)7^5^ 
2/cp 


+00 


- Pi - jfcp)fcp, 


( 6 ) 


]=-oo 


where the summation goes over all modes with frequencies pj = Pi + [j — l)/cp, and is extended 
to ±00 since the total width of the rf excitation comb APcomb is much larger than /cp and the 
homogeneous linewidth Appom. The {I — m){I + m + 1) factor arises from the dipolar transition 
matrix element!^. We note that Supplementary Eq. [^does not involve any explicit nuclear-nuclear 
interactions. Instead such interactions are introduced in Supplementary Eq. [^phenomenologically 
via the homogeneous broadening described by the lineshape function L{v). On the other hand the 
presence of finite homogeneous broadening is essential in order to use the limit of weak rf field^^ 
and transform the Bloch equations into rate equations (Supplementary Eq. j^. The validity of the 
weak rf field approximation is discussed and verified experimentally in [Supplementary Note 3 


Since Supplementary Eq. is a system of linear first-order equations, the solution is a multi¬ 
exponential relaxation towards the fully depolarised state where all nuclear spin states have equal 
populations pm = 1/(2/ -|- 1). The solution has the general form: 

21 

Pm{t) = 1/(2/ + 1) + ^ CLm,] exp(-Ajt), (7) 

j=l 


where Aj are the non-zero eigenvalues of the ODE system matrix of Supplementary Eq. [^ The 
values of Aj depend on all transition rates Wm,m+i from Supplementary Eq. [^ while the coefficients 
am,] depend both on Wm,m+i and the initial probabilities pm,o from Supplementary Eq. [^ 

Non-zero nuclear spin polarisation along the magnetic field {Oz axis) changes the spectral split¬ 
ting of the quantum dot Zeeman doublet. Such change known as the Overhauser shift is measured 
experimentally using photoluminescence spectroscopy. The time evolution of the Overhauser shift 
for the fixed values of nuclear transition frequencies {Pm,m+i} = ■■■ 

reads as: 


+/ 

/'hf,!(i)/cP))/nuc)/'(r^)) — -d ^ ^ (8) 

m=—I 

where A is the hyperfine constant and we used Supplementary Eqns. Bin so that Ehf,! is 
dependent on /cp, Bi, T/uo the homogeneous lineshape function L(p) and all nuclear transition 
frequencies {i2m,m+i} as parameters. 

Each QD contains a large number of nuclear spins with randomly distributed absorption fre¬ 
quencies. Thus to describe the experiment on nuclear spins in a self-assembled quantum dot we 
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need to average over all 1 ^ 111 , 10 + 1 , which can be done over one period /cp since the spectrum of the 
frequency-comb rf excitation is periodic. Similarly to the case of / = 1/2, the following expression 
is obtained for the time dependence of the Overhauser shift, describing the dynamics of rf-induced 
nuclear spin depolarisation: 


Ehf{t, fcP, Bi,Tnuc, L{i')) = 

(9) 

f(jp I ... I duj—i j iHpf 1 (t,/cp, Hi, T/oq, Z/(p), {Poi,m+l})• 

0 0 

The quantity measured in the experiment is the rf-induced variation of the Overhauser shift: 


/cP, -Bi, Tnuc, T(p)) = Ehf{t, f CP, Bi,Triuc, L{iy)) — £'hf(t = 0). 


( 10 ) 


The values of fcp and Bi are the parameters that are controlled in the experiment. The nuclear 
spin temperature Tnuc can be determined using the known hyperfine constant A and the measured 
total Overhauser shift A£'hf(t = 00 ). Thus for a given homogeneous NMR lineshape L{v) the 
nuclear spin depolarisation dynamics can be fully predicted from Supplementary Eqns. [3' 10 Con¬ 
versely, Supplementary Eq. [^can be treated as an integral equation on the unknown homogeneous 
lineshape function L{v). Since Fredholm’s integral equation of the first kind is an ill-conditioned 
problem, some constrains on L{v) are required. Our approach is to use a model lineshape with 
two parameters Aphom and k (Eq. 1 of the main text). Upon substituting this model lineshape 
the Overhauser shift variation of Supplementary Eq. [^becomes AEhf(t, /cp, .Bi, Tnuc, Aphom, k). 
We then perform least-square fitting to the experimental dependence A£'hf(t,/cp) using Bi, the 
homogeneous linewidth Aphom and the roll-off parameter k as fitting parameters and using the 
nuclear spin temperature Tnuc obtained from the experiment. 

The ODE system of Supplementary Eq. can be solved analytically for I < 3/2, however 
it turns out to be more practical to perform numerical diagonalization in order to obtain the 
eigenvalues Aj and coefficients am,] which are then used in Supplementary Eq. Similarly we use 
numerical integration to evaluate Supplementary Eq. 


Supplementary Note 3. APPLICABILITY OF THE FREQUENCY COMB TECHNIQUE 
AND THE RATE EQUATION MODEL. 

The evolution of the nuclear magnetization under rf excitation can be described by the Bloch 
equation^^. In this model the solution under resonant monochromatic excitation is determined by 
the three important parameters: the amplitude of the resonant field Brf, and the relaxation times 
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characterizing the system, the longitudinal Ti and the transverse T 2 . In self-assembled quantum 
dots Ti is extremely long (few hour^ go that the longitudinal relaxation can be neglected. 
Thus the nuclear spin dynamics is determined by the relation between Bj-f and T 2 . Two cases are 
possibld^. If rf magnetic field is strong {'yBr{T 2 S> 1) the nuclear magnetization has oscillatory 
behaviour (Rabi oscillations are observed). By contrast, for weak rf excitation {jB^{T 2 <C 1) there 
are no oscillations, and any nuclear magnetization along the external field decays exponentially 
to its steady state valu^^. The exponential dynamics in the weak rf excitation regime allow for the 
problem to be simplified and for the rate equation model described by Supplementary Eqns. [^to be 
used. The validity of the rate equation model is essential for the determination of the homogeneous 
lineshape and thus sets the applicability limit for the frequency comb technique itself. 

To verify the applicability of the rate equation model we performed frequency comb spectroscopy 
measurements at different rf amplitudes. Supplementary figure shows the results for ^^Ga mea¬ 
sured at high rf field density /3i = 86.0 nT Hz“^/^ and low rf field density /3i = 24.2 nT Hz“^/^ - 
in all other respects the conditions in these experiments were the same as in the experiment with 
medium /3i = 39.1 nT Hz“^/^ shown in Fig. 2 of the main text. 

The Supplementary figure shows by symbols the nuclear spin depolarisation dynamics at 
small comb period /cp = 31.25 Hz. At low rf field amplitude /3i = 24.2 nT Hz“^/^ (triangles) the 
decay can be described very well by a single exponential decay (r ~ 290 ms) shown with a solid 
line. By contrast, at high amplitude j3i = 86.0 nT Hz“^/^ (squares) the dynamics shows signatures 
of oscillations, and there is a clear deviation from the exponential behaviour (best exponential fit 
is for r ~ 23 ms). 

Supplementary Figures and c show the full AFhf(t, /cp) dependencies measured at high 
(b) and low rf (c) field densities. The AFhf(t,/cp) profiles are in good agreement for the two 
experiments (except for the rescaling along the trf axis). Using model htting we find Aphom ~ 
244 Hz, k « 1.66 for f3i = 86.0 nT Hz“^/^ and Ai^hom ~ 223 Hz, k « 1.68 for /3i = 24.2 nT Hz“^/^. 
This is also in good agreement with Aphom ~ 221 Hz, k ~ 1.67 found for the measurement at 
/?! = 39.1 nT Hz“^/^ shown in Fig. 2 of the main text. 

We thus conclude that the discrepancy between the results in Supplementary Figs. and c 
becomes significant only at small comb period /cp = 31.25 Hz (as also demonstrated in Supple¬ 
mentary Figure]^). This can be explained as follows: At high rf excitation amplitude the nuclear 
spin depolarisation takes place on a shorter time scale r. If the rf pulses are shorter than l//cp, 
the spectral profile of the frequency comb becomes distorted. Thus if the nuclear polarisation 
decay timescales r are shorter than l//cp, the rate equation model is no longer applicable since 
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the rf excitation can not be described as a frequency comb. Thus it is required that r > l//cp. 
Furthermore, to measure the homogeneous lineshape and linewidth APhom we only need to use 
frequency combs with comb periods /cp comparable to or larger than Aj/hom, so it is required that 
fcp G {APhom) oo}. Combining r > l//cp and /cp G {APhomjOo} we find the following condition 
on the frequency comb technique applicability: 

APhom > l/r, (11) 

which restricts the rf amplitude, characterized by the depolarisation time r. 

Another requirement, arising from the applicability of the weak rf field limit of the Bloch equa¬ 
tions, is that the rf induced depolarisation time r must be longer than the transverse relaxation 
time T 2 . However, T 2 is related to the homogeneous linewidth as T 2 ~ l/(7rAPhom)- Thus the 
requirement r > T 2 leads to the same condition as that of Supplementary Eq. Furthermore, 
since the frequency comb technique relies on the measurement of the longitudinal nuclear magne¬ 
tization, the depolarisation time r must be shorter than the nuclear spin Ti times. In combination 
with Supplementary Eq. [m this leads to the following condition: 

APhom > 1 /t > l/7"i (12) 

This condition has a dual role: it sets the boundaries for the rf excitation amplitude (characterized 
by r) and sets the limitation Aphom > 1/T"i on the properties of the nuclear spin system that can 
be studied with the frequency comb technique. This latter condition can be rewritten as 

Ti > aT2 (13) 

From the measurements at different rf amplitudes we find that the parameter a must be a ~ 10—100 
or larger in order for the frequency comb technique to work reliably. This however is a rather weak 
condition and is satisfied for a large class of solid-state nuclear spin systems where Ti ^ T 2 . This 
demonstrates the wide applicability of the frequency comb spectroscopy technique developed here. 
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